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Abstract. The so-called "phenomenological" kinetic equation for one-pair density operator 
for spin-selective reactions is defended. We derive this equation from the kinetic equation for 
density operator of all pairs which are treated as singlet and triplet bosons. There presented 
£^ ' some reasons for inconsistency of measurement-like approach to the problem. 
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C«V 1. Introduction 

On ' 

We address the problem of correct treating and interpretation of recombination events in the 
ensemble of radical pairs. Particularly we deal with the relation between the spin-selective 
recombinations of pairs and measurements of their spin. Our main statement concerns the 
radical distinction between a measurement-like process and proper act of recombination. In 
numerical works by Kominis (see [1J and references therein) any recombination event is 
considered as that of measurement-like. On the basis of this view-point his critics of the 
well-known 'phenomenological' kinetic equation JHO 

jP + i\H,p] = -k s (Q s p + pQs) ~ k T {QTp + pQ T ) (1) 

is based (h = 1). Here p is the density matrix of a pair; H - its Hamiltonian; 2k$ and 2kx are 
the rate constants of singlet and triplet recombinations respectively. For simplicity we take 
into account only spin variables of the pairs. In this case p is an operator in a 4-dimensional 
Hilbert space with the basis {\S), |T_), \T ), |T+)}; 

Q s = \s)(s\, q t = J2 l T -)^l < 2) 

<T=0,±1 

are the projectors onto the corresponding spin subspaces. In literature one can find a 
misleading interpretation of the last two anticommutator terms in rhs of (0Q) as remnants of 
a full Lindblad-type structure [4] upon elimination of the so-called 'sandwich' gain term: 

2QpQ-Qp- P Q ^ ~(Qp + pQ) (3) 
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(here Q stands for Q s as well as for Q T ). In accordance with this interpretation the 
transformation © is considered as the elimination of the gain term because of recombination. 
Note that the Lindblad-type structure in lhs of © provides strait associations with the spin- 
measurement-like process. The central point of the Kominis' approach is the restoration of 
the full Lindblad structures in © in the framework of the measurement paradigm. 

We are going to show that there is a firm background behind the equation ©, which has 
practically nothing common with any measurement-like process. 

2. Spin-selective recombinations 

Our starting point is the kinetic equation for the multi-pair quantum state g. In 
contrast with p this operates in the infinite-dimensional Hilbert space with the Fock basis 
{\ns, riT + , tit , where n are the numbers of pairs in the corresponding states. We 

should write the kinetic equation in terms of annihilation, as, a TtT , and creation, a' s , Oj, , 
operators for singlet and triplet pair respectively. So we treat any pair as an individual 
'particle'. These 'particles' are bosons, i.e. one has the following commutator relations: 

[«5,4] = !; [flr ri .4j = <W2! [05, 4j = °- ( 4 ) 
The mathematical description of the process of spontaneous acts of death of the 'particles' (the 
acts of recombinations) is known. For the both types of recombinations (singlet and triplet 
ones) we have the corresponding three-term Lindblad structure in rhs of the kinetic equation: 

-^p(t) + i[H, g(t)] = k s (2a s g(t) at, - a< s asg{t) - g{t) a^a 5 ) + (5) 
k T ^2(2a Ta g(t) a ] Ta - a ] Ta a Ta g{t) - g(t) a^a T J 

The Hamiltonian H in © should not be confused with H from ©. The evaluation of this 
equation is made in Appendix. 

The one-pair density matrix p from © is obtained by the following standard way: 

p = \S)(ala s )(S\ + (\T„)(ala T J(S\ + \S}(ala s )(T a \) + (6) 

l T *l)( a T CT2 «T CTl )(r<x2|, 

0-1,0-2 

where (...) = Tr(g(t)...). The mean value of any multi-pair operator O evolves according to 
the equation 

j t (6) + i([6,H}) = k s ((al[6,a s }) + ([ald}a s )) + (7) 

^E(( a T CT [o,a T j) + ([4 CT ,d]a TCT )) 

a 

Substituting the operator O with proper bilinear products of a) and a, one gets the kinetic 
equation for all elements of p: 

^-(a ] s a s ) + i({ala s , H]) = -2k s (ala s ) (8) 
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4(4«T ff ) + i{[a s a Ta ,H\) = ~{k s + k T )(a s a Ta ) (9) 
^(«t ct «s) +i{[a i TiT as,H}) = -(k s + k T ) {a ] Ta a s ) (10) 

^(«L«t ct1 ) + «([a^ 2 a TCTl , H\) = -2k T {a ] Ta ,a Tal ) (11) 

Combining these equations with the use of © and © we arrive at an equation for p with rhs 
identical to that of ® • 

Its worth to stress that the equation © describes the process of spontaneous spin- 
selective 'annihilation' of pairs, not measurement of their spins. 'If it were the case' picture 
of a measurement-like process will be discussed later. 



3. Hypothetical measurement-like recombinations 

Let us consider an imaginary process in which recombination and pair's spin measurement 
are really inseparable. It will be shown that the evolution equation for the one-pair density 
matrix p resembles neither (OQ) nor the equation by Kominis HI. 

Note that any measurement of a pair's spin create information about the spin value. This 
information is written down in the environment. Hence the pair appears to be entangled with 
the environment. Let us first consider the following phenomenological equation for p where 
the pairs are subjected to spin measurements but no recombinations take place: 

j t P + i[n,p] = k(Q sP Q s + Q T pQ T -p). (12) 

Here k is the number of measurement events per second. The first two terms in the brackets in 
rhs of (fT2l) stand for the two possible outcome of a measurement (singlet or triplet states). Any 
singlet-triplet coherence vanish in accordance with rhs of (OQ). This is due to the mentioned 
entangling process. At the same time the total number of pairs conserves: Trp = const (no 
recombinations). Now assume that upon a spin measurement with the singlet outcome the 
pair recombines (in the singlet channel) with the probability ps, and the same for the triplet 
outcome. One may account the loss of pairs due to these recombinations eliminating parts of 
the gain terms in rhs of (fT2)) : 

j t p + i[U, p] = fc((l - Ps)QsP Qs + (l~ Pt)QtP Qt-p). (13) 

If we introduce the notations k s = psk, kr = Prk and take into account that Q T = 1 — Qs, 
we arrive at 

j t P + i[U, p] = (2k -k s - k T )QsP Qs~{k- k T )(Q sP + pQs)~ k T p.(U) 

The 'sandwich' as well as anticommutator terms are present. Note that if k T = and k s = k 
the equation (fT3l) lies just in between the same cases of (OQ) and equation (5) from (TJ. 

In the work [5] the equation (fl4)) for the special case k = ks + kx was developed. The 
comparative numerical analysis of the derived equation and that one from (OQ) revealed a slight 
difference. 
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4. Conclusions 

Resuming we state that the kinetic equation for one-pair density matrix p in the case of spin- 
selective recombinations has been shown to have the following form: 

= -i[H,p] - k s {Q s p + pQs) - k T (Qrp + pQt)- 

Our approach is only valid for geminate pairs. In the opposite case (for possibility of cross- 
recombinations) one may not consider the radical pairs as bosons. 

We show that the kinetic equation for the measurement-like recombination process 
differs radically both from the conventional equation and from Kominis' results. We also do 
not share the opinion from [[51 "...that the the quantum measurement approach should normally 
be used in future simulations of spin-selective radical pair reactions". From our view-point the 
spin selective recombination is not a measurement. In the hypothetic opposite case when any 
recombination is preceded by spin measurement there emerges an information concerning the 
spin of survived pairs. This options corresponds to the terms from (fT3l) proportional to 1 — ps 
and 1 — p T . There is no any basis to suppose such a phenomenon takes place in recombination 
process. 



Appendix 

There is a regular derivation of equations like ([5]) from the first principles j6l 13- As an 
example, we will give a sketch of such a derivation for singlet Lindblad structure in rhs of 
©. For example it can be obtained in the second order of perturbation with respect to the 
following dynamics: 

V = \a s ® (b f -E 1 ") +H.c. (15) 

responsible for transformation of a singlet pair into a product b. This process is accompanied 
by photon emission. Electromagnetic field is the vector operator E* (this part is responsible 
for photon creation). So the product should be vector as well, to make V invariant under 
rotations, i.e. it should have total unit angular momentum. Let R(t) be the total statistical 
operators of pairs, products and photons. R(t) evolves with respect to the following equation 

jR{t) = -i[H a + H b + H ph + V, R(t)}, (16) 

where H a , H b and H p h are free Hamiltonians of pairs, products and photons, respectively. In 
the interaction representation 

R^i) = exp[i(H a + H b + H ph )t)R{t) exp[-i(H a + H b + H ph )t], (17) 

Vj(t) = exp[i(H a + H b + H ph )t]Vexp[-i(H a + H b + H ph )t] (18) 

(TT6b turns into 

^R I (t) = -i[Vi(t),R I (t)]. (19) 
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In the second order with respect to Vi(t) one has 

/t+At 
tft'tVXO, - (20) 



t+At 



df / dt'ivtf), [Vj(nRi(t)}} 



In the last term in rhs we replaced Ri{t") with Ri(t) due to the assumption At <C r re/ , where 
r re i is a typical relaxation time. There is another time scale, r corr - the correlation time which 
will appear later. If At r corr one can neglect all correlations between pairs and their 
environment (products and photons): 

Rl(t) = Qi(t) ® Qb,ph (21) 
The statistical operator g btPh is assumed stationary. We also suppose that 

Eg b , ph = Qb, P h& = 0, (22) 

i.e. all emitted photons are rapidly absorbed so that g biPh is the vacuum with respect to 
electromagnetic field. 

From (l20l) we get the 'coarse-grained' equation 

AQ!(t) _ 1 rt ' 

At 

Omitting details we arrive at 

\I2 poo pt+At 

g{r)(aUtf)a s {l/-T)g I {t)- (24) 



i pt+At pt' 

-— J dt' J df'Tr^V^t'), \Vi(f),gi{t) ® g b , ph }}. (23) 



A m |Ap r [»* 



At At J 

a s (t' - r) gi (t)al(t')^ + g{-r) (^(t)a^(t' - r)a 5 (t') - a s (t')M*)4(*' " r ) 
Here 

ff(r) = ff(-r)* = Tr b>ph (bi(r) • Ej(r))(b{(0) • E{(0))^ (25) 

is the correlation function; r corr is its typical scale. For simplicity we assume (H a )r corr <C 1. 
Coming back to the Schrodinger representation we get 

lg(t)+i[H a ,Q(t))= (26) 

(k s + iK S )(a s g(t) al - a s asg(t)) + (k s - iKs)(a s g(t) a s - g(t) agas), 

where 

POO 

k s = \X\ 2 Re / g(r)dT (27) 
Jo 

is the rate constant of singlet recombinations and 

POO 

k s = \\\ 2 Im / g(r)dr (28) 
Jo 

is a parameter of a slight Hamiltonian renormalization: 

H a ^H a + Ksa j s a s = H. (29) 
We get the first line of ©. 
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